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Abstract 



We review the relation between the matrix model and Liouville approaches to two- 
dimensional gravity as elaborated by Moore, Seiberg and Staudacher. Then, based 
on the supersymmetric Liouville formulation and the discrete eigenvalue model pro- 
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tion to the supersymmetric case. The minisuperspace approximation for the super- 
symmetric case is formulated, and the corresponding wave equation is found. 
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1. Introduction 

It is generally accepted that the m-th multicritical point of the Hermitian 1-matrix 
model corresponds to the coupling of the Liouville theory to perturbations of the 
(2, 2m — 1) minimal conformal model. However it took some time to understand 
this relation. A precise dictionary between Liouville and matrix model operators, 
correlation functions and loop amplitudes was proposed by Moore, Seiberg and 
Staudacher in Ref. [1]. They recognized how part of the difficulties were related 
to contact terms. After a careful examination of macroscopic loop amplitudes in 
two-dimensional gravity, which must satisfy the Wheeler-de Witt (WdW) equation, 
they managed to determine a new frame of scaling operators and couplings (called 
conformal frame) in the matrix theory, whose correlation functions and scaling 
dimensions were in perfect agreement with the Liouville predictions. 

This successful program motivated us to generalize their computations for the 
super symmetric case. We considered the amplitudes of the N = 1 super-Liouville 
theory coupled to c < 3/2 matter, which were calculated in Ref. [2]. As for 
the discrete counterpart, there is no actual generalized matrix model. So far the 
most successful description has been given by the discrete model introduced by 
Alvarez-Gaume, Itoyama, Mahcs and Zadra in Ref. [3] in terms of super-eigenvalue 
variables. Nevertheless we shall prove that, in spite of the missing interpretation 
in terms of matrices, this discrete theory also shows a complete agreement with 
continuum results. 

We dedicate sect. 2 to a brief review of the purely bosonic case, emphasizing 
the main results from the continuum approach in subsect. 2.1 and describing in 
subsect. 2.2 the discrete theory and the procedure that determines the conformal 
frame of operators. In sect. 3 we present the super symmetric version of the models: 
the super Liouville theory is summarized in subsect. 3.1; the super-eigenvalue 
model, its standard scaling operators and their correlation functions are reviewed 
in subsect. 3.2; then, in subsect. 3.3 we follow the strategies of the bosonic case 
to define a superconformal frame; their wave functions are studied in subsect. 3.4. 
We leave sect. 4 for the final comments and a brief conclusion. 

2. Review of the purely bosonic case 
2.1 Continuum approach 

Let us take the David-Distler-Kawai^ (DDK) description of c < 1 conformal theories 
coupled to 2d gravity, where the partition function Z on the Euclidian sphere is 
defined by the following integral over "matter" {X) and Liouville ((/>) fields, 

^=/™e-, (1) 

= ^ y e?w^[g^''da^di,^ + r%XdbX - QR(t> + 2iaoRX + 87r/2e"'^] . (2) 

Qab as in (2) is the fiducial metric fixed by the conformal gauge, R is the respective 
scalar curvature and Vsl{2,c) is the volume of the residual SL{2, C) symmetry (a 
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residual symmetry of the conformal gauge). The parameter ao and the background 
charge Q are related to the matter central charge c by the relations 



l-12a^ , Q = 2^2 + al^^'^^ . (3) 

The last term in Eq. (2), // J d^u; \fgc°"f', measures the area of the spherical surface 
and is called the cosmological term. The coupling /i is a bare cosmological constant 
and the parameter a is chosen in such a way that e""^ has conformal dimension one, 

A(e«^) + g) = l . (4) 

The above equation has two solutions, namely a± = —Q/2 ± |q;o|. We choose 
a = a_i_ > —(5/2, which exhibits the correct classical limit and also satisfies Seiberg's 
energy bound^ {E{a) = a + Q /2 > 0) . 

As stated in the Introduction, we are interested in the scaling behaviour of 
correlation functions of vertex operators. In the continuum approach we take vertex 
operators "dressed" by gravity, Vi = f \fgO iC^^^"^ ^ where Oi is the matter field 
or vertex and the dressing factor e^*"*^ is determined by imposing that OjC^*'^ has 
conformal dimension 1, 

A(C»ie^*'^) = /ii-^/3i(A + g) = l . (5) 

Above, hi is the conformal dimension of Oi. In particular, for a tachyon Oi = e*^*^, 
one finds hi = ^ki{ki — 2aQ) and therefore its dressing charge reads (3i{ki) = —Q/2+ 
\ki — ao\. Notice the degeneracy in this representation: the momenta ki and 2aQ — ki 
imply the same conformal dimension and dressing. Therefore we may take either 



(6) 



to represent the same dressed operator, choosing one or the other according to kine- 
matic criteria. We call the above vertex (T^) the tachyon, and k is its momentum. 
In particular the cosmological term or area operator Tq = f d^z y/ge°'+'^ previously 
mentioned may be as well represented by T2ao— I Vq e^*'*°"'^"'"'*+'^. 

An A'"-tachyon correlation function An = {TkT_ ■ ■ - Tk^^) is calculated, using 
Coulomb gas techniques. We choose the fiducial metric ga/3 = Sai3- Integrating over 
the matter field X we obtain the momentum-conservation law 

N 

X!^i = 2«o , (7) 

which must be satisfied by every non-vanishing correlation function. 

The gravitational contribution is more difficult to deal with, due to the ex- 
ponential self-interaction in the Liouville action. At this point we turn to the 
zero-mode integration procedure^: defining (f) = cj) + (f)Q, where 0o is the Liouville 
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zero mode, we can explicitly integrate over (I)q and use the free field propagator 

(e/^0(^)e7<^(^)) = {z- w)-f^^ for the remaining field 4>. 

After renormalizing the cosmological constant {p, —>■ ji), we find the A'"-point 
function 

where A(x) = and, = -^(Q + E.li ft) ■ 

After writing Eq. (8) we analytically continue the results, originally derived 
for s e 2, to s G (D. Moreover we constrain the kinematics choosing N — 1 values of 
ki to be bigger than and one smaller (for a thorough discussion on the subject 
see Refs. [7,10]). Notice that An factorizes in terms of external leg contributions 
and its scaling behaviour with respect to the cosmological constant ~ = 

jj^-Q/<^+-J2i=i ^ provides the scaling dimension of the dressed vertex operators 

Tfc = y'd2^e^'=^+^(")'^ ^ . (9) 

The scaling of the partition function can also be taken from Eq. (8) and defines 
the so-called string susceptibility 7 as 

Z ~ A ~ //"^/"+ = , (10) 

i.e. 7 = 2 + Q/a+. 

From the momentum conservation (7) we conclude that there is only one tachy- 
onic operator, with momentum k — 2ao and dressing charge /3 = a+ (i.e. the area 
operator), whose expectation value is non- vanishing, 

(T^ao) ~ . (11) 

The remaining 1-tachyon functions are expected to vanish. 

The conservation law (7) also implies that T^^ and T2ao-ki represent the 
same operator, and we conclude that 2-point functions are diagonal and scale as 

^-Q/Q+-2/3(fc)/a+_ 

In Ref. [1] the authors used the wave functions of dressed vertex operators to 
establish the correct Liouville/matrix model dictionary. The starting point was the 
WdW equation^ 

[H-Ao]^i{O) = , (12) 

where the wave function \I/;((9) corresponds to the insertion of a dressed operator 
C on a surface of boundary I; Aq = 1 — h is the conformal weight of C; H is the 
Liouville Hamiltonian 



= l(,6' + 47rP)2 + Q(0' + 47rP)'+-^e«+^+^ . (13a) 
2 2q!_i_ 8 
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In the minisuperspace approximation one ignores the space dependence and drops 
derivative terms in which then reads 



with p = —id/d(j). Moreover, since the boundary length I is measured by the contour 
integral Z = ^d^ e2"+'^, we can associate 



d<P 2\!:dl) ' 2al^ 8 ^' ^ ' 



in which case Eq. (12) becomes the Bessel differential equation 



*KO) = , (15) 



where = = ( ^ + . 

Since the wave functions '^liO) are expected to decay in the infrared (large V) 
limit, we take the modified Bessel functions of second kind \I/z(C) oc K^{y^l) as 
appropriate solutions. This result will be taken as a guide to unravel some tangles 
in the continuum/discrete translation. 

2.2 Discrete (matrix model) approach 

The Hermitian 1-matrix model is defined by the partition function 

Z = j V(j)e-^''^^'t'^ , ^ (16) 

i i<j k 

One can generate different critical regimes by tuning the coupling constants gk- The 
potentials found by Kazakov in Ref. [8] exemplify this property. In the m-th critical 
regime the continuum limit of the model is defined^ by a double scaling limit, where 
the matrix size gets large {N — > oo) and the constant A approaches a critical value 
(A Ac), while the combination A?" (Ac — A)^'^^/^'^ is kept constant. The resulting 
theory is described by the free energy 

^=-^di^Mtn) , (17) 

where k is the renormalized string coupling and u is the specific heat satisfying the 
string equation 

-to = '^ntnRn[u] ; (18) 

n=l 
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i?„[it] arc the Gel'fand-Dikii polynomials. The dependence on the renormalized 
couplings tn is ruled by the KdV flows 

d d 

u=-Rn+i[u] . (19) 



dtn dt 

The exact m-th critical regime is defined by the limit 

to^t ; ; ^^^^ 

tn ^ Q ) n ^ 0, m , 

in which case the string equation becomes mRm[u] = t. Within the planar or 
spherical approximation, given by the limit k ^ 0, the Gel'fand-Dikii polynomials 
tend to Rn[u] — u^/n, and the string equation (18) becomes 



J]^n«" = , (21) 



n>0 

which in the m-th critical regime (20) implies = t. 

We associate a set of scaling operators (t„ to the couplings i.e. cr^ ^ 
This frame of operators and couplings was called the KdV frame by the authors in 
Ref. [1]. Correlation functions are described by the KdV flows, with the following 
results in the planar approximation 

n 




, a = J2^i ■ (22) 

In the regime (20) we observe the scaling behaviour of the free energy as given by 
JF ^ f'2+i/m _ ^2-7^ which defines the susceptibility 7^, = — 1/m. On the other 
hand the n-point functions scale as ^2+i/m+^(ai/m-i)^ which therefore defines the 
scaling dimension dn — n/m oi the operator an- Notice that neither any one of the 
1- nor of the 2-point functions vanish, as opposed to the continuum results. Such 
differences obstruct a direct association between vertex operators and the KdV 
frame of scaling operators an- 

These difficulties can also be felt as one studies the wave functions of scaling 
operators. We start defining macroscopic loop operators W{1) of finite length Z, 
whose expectation value reads 

{W{1)) = —= \ dy e-''^^^'*") . (23) 

We suggest Ref. [8] for a review on the subject. Notice that we can actually 
integrate-*^ Eq. (23) by using the planar string equation (21), finding as a result 

(^(;)) = J_ Vt„«'^+^/Vn(M , (24) 
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v., 



with 

n{x) = ^ dz (-) e-^ = n\x—'/'e-^ J] ^ . (25) 

The functions i/^ni^u) display the dependence on the loop length I. They are pro- 
portional to the wave functions associated to the insertion of scaling operators, as 
follows 

*/(^n) = {(TnW{l)) = -^{Wil)) « Vn(/W) , (26) 

i.e. Eq. (24) can be taken as an expansion of the loops in terms of wave functions. 
However we see that il^nilu) does not obey a Bessel equation, which reinforces the 
differences between the KdV basis and the vertex operators. 

In Ref. [1] the authors resolved these discrepancies defining a new frame of 
scaling operators dn and the respective couplings tn, which they called the confor- 
nial field theory frame. They assumed that the differences between the observed 
correlation functions should arise from contact terms. On the other hand a change 
of contact terms is equivalent to an analytic redefinition of coupling constants. Thus 
one should look for the proper set of couplings. 

In order to change coupling constants and yet preserve scaling properties we 
must examine their dimensions. From the way the loop length / comes out in 
Eq. (23), we take the dimension of the specific heat, [u] = [length]"^. Thus the 
dimension of the couplings follow from the string equation (21), 

[y = [length]— . (27) 
The "physical" cosmological constant is usually taken to be the one coupled to 
the area operator and should therefore have the dimension of inverse of area or 
[length]"^. Since [t] — [to] — [length]"'" we conclude that t deserves the name 
of cosmological constant only when m = 2. In general tm-2 is the coupling with 
the proper dimension of a cosmological constant. If one wishes to compare results 
with the Liouville approach one should select another regime, different from the 
one defined by (20), where scaling properties are measured with respect to some 
coupling as tm-2- 

We can determine the conformal frame of couplings according to the following 
criteria. Consider an analytical transformation relating the KdV frame {tn} to 
another frame of couplings {tn}, with the general form 

tn = ClU + Cii%% + • • • (28) 

In order to maintain scaling properties (the operators dn should be scaling fields 
as cr„ already were), let us impose that the analytic transformation preserves the 
dimension of couplings, [t^] = \tn\- In analogy to the limit (20) we shall define 
another m-th critical regime in terms of the new set of couplings, 

tn ^ , Uy^ 171 — 2,171 , 
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where /i is the "physical" cosmological constant. The solution of the string equation 
will become v? = ^ and correlation functions will scale as functions of The 
operator am-2 is expected to correspond to the area or cosmological term. 

Our immediate concern is to have 1- and 2-point functions compatible with 
the Liouville predictions. In this case it is sufficient^ to consider lowest (linear)- 
order terms in the transformation (28) around the critical point (29). Thus it is 
convenient to shift the non- vanishing couplings as tm-2 —>■ A* + ^m-2, tm —>■ — 1 + ^m, 
so that the critical regime is achieved by the limit tn ^ and we can deal with the 
tnS as perturbative couplings. 

Considering these scaling restrictions and lowest-order approximations, the 
transformation (28) is reduced to the form 

oo 

tn = + J2 oi"+'^^A*'tn+2. , (30) 

where 6„ and a^J^^ are dimensionless coefficients to be determined. However, the 
coefficients 6„ are not independent: they are given by 

bm-2i = tO-!?i~^'' ; bn^O , n> m . (31) 

From Eq. (30) we derive 7^ = Yl^I^ln^ af^^ jrr^ — , which gives the transformation 

Otn — OZrL — 2s 

of scaling operators 

[n/2] 

an=Yl «^^«''^n-2s . (32) 
s=0 

(n) 

Therefore the coefficients as uniquely characterize a given basis. In order to find 
the conformal basis we can adopt two strategies, described below. 

2. 2. a. Minisuperspace approximation strategy 

In this case we use the wave functions as a guide-^ to determine the coefficients a^J^\ 
From the Liouville description (see Eq. (15)) we expect to find Bessel functions. 
We take the expansion (24) for the 1-loop function, which holds in any regime, and 
substitute the transformation (30). Recalling that — > /x we find 

m/2] 



{W{1)) = -^U^+'/^ V bm-2si^m-2s{lu) 



s=0 

/2] 



(33) 



+ ^ J]tn«"+'/' E «i"^V'n-2.(i«) 



n>0 \ s=0 

Indeed we can obtain Bessel functions using 

^/2] 



Kn+l/2{x) = a / 2 5Z a^J'^fPn-2s{x) 



(34) 



7 



with 

' ~ 2i s\{i-s)\{i-2s)\ ■ ^ ^ 

These are the coefficients we were looking for. Using Eq. (31) we also find the 
corresponding expression for the remaining coefficients 

[m — 1/2) (m — 1/2) 

Therefore the 1-loop function becomes an expansion in Bessel functions 

^^^^^^ = ^ {m-112) \-^rn-V2{lu) - i^n^+l/2(H] 

x/2 ^37) 

+ ^ Vrn«"+^/'i^n+l/2(M , 
KIT ^ — ' 
n>0 

from which we obtain the wave functions 

*K^n) = {^nW{l)) = ^u^+^/^K^+y^{lu) , (38) 

KTT 

which obey the Bessel equation (15) for u = n + 1/2. Notice also that in the exact 
regime (29) we have 

./o ym+1/2 

= -y— «™-'/'i^m-i/2(M = -y(?m-iW-(0) , 

i.e. (am_iVF(/)) = — |(VF(Z)), indicating-*^ that ct^-i is the boundary operator 
which measures the loop length I. 

We remark that the coefficients (35) of the conformal basis have an interesting 
interpretation in terms of orthogonal polynomials: indeed from Eqs. (25) and (34) 
we can write 

where Pn{x) = X^ilfo^ ai"'*x'^~^* are the Legendre polynomials. We shall understand 
the appearance of this family of polynomials in the following subsection. 
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2.2.b OHhogonalization strategy 

(n) 

The coefficients (is can also be determined by imposing that the 2-point functions 
be diagonal. In fact, from the general transformation (32) and the 2-point 
functions of the KdV frame (aiaj), which follow from Eq. (14), in the regime 
ji, we find 

^(i+j+l)/2 

{aiaj) = ^2 , (41) 

where 

[i/2][j72] (i) (i) 

can be seen as a characteristic metric of the basis defined by some coefficients ai"^^ . 

Our aim is therefore to choose ai"^-* so that gij becomes diagonal. Observe that gij 
satisfies 

[1 _ (-1)^+^+1]^. .(a) = / dxPi{x) Pj{x) , Pi{x) = J2 4^^'-^' > (43) 

•^-1 s=0 

where Pi{x) = Yl^s=o are characteristic polynomials of a given basis. Thus 

our problem has been translated into finding, among the characteristic polynomi- 
als, those which are orthogonal with respect to the internal product {Pi, Pj) = 

J^^ dx Pi{x)Pj(x). The well-known solution is given by the Legendre polynomi- 
als, whose coefficients are precisely those defined in (35). Therefore the 2-point 
functions in the conformal frame read 

= j-i^rr^^j- ' i + ieven ^44) 

t analytical in , Z + J odd 

i.e. they are diagonal up to analytical terms in the coupling constants, as expected 
from the Liouville approach. 

Following any of the strategies we conclude that the conformal frame of oper- 
ators an is related to the KdV frame an by 

^ [^l(-l)^ {2n-2s)\ , ,^ , 

2^ s\[n — sy\n — 2s)\ 



which can also be inverted 

[n/2] ^ 

"2" s!r(n-s + 3/2) 



-r;^ v^n! (2n-4s + l) 



s=0 
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We must also test the 1-point functions. From Eqs. (44) and (456) we can 
calculate ^(cn) = which, after integration, gives 



{a. 



m-2+2s 



) = 



~^ 2m-2+2k k\T{m + k + 1/2) 



Now we can take Eq. (45a) and easily calculate the 1-point functions 
out that the only non- vanishing cases are 



-1/2 



Av2 (2m- 1) (2m -3) 



C^m) = 



m+1/2 



k2 (2m + l)(2m - 1) 



(46) 



It turns 



(47a, b) 



This is also in agreement with the Liouville approach: {am-2) is the average area, 
as in Eq. (11); on the other hand (am) can be interpreted^ as a non-tachyonic 
part of the average energy. Since {am-2) = ^-^7 we can integrate Eq. (47a) 

and obtain the scaling behaviour of the free energy in this regime: J-' ~ fi^^^^"^. 
Higher-order correlation functions would require the analysis of higher-order terms 
in the transformation (30). Nevertheless we know that, within the regime (29), the 
A'"-point functions must scale as 



N 

n 



m+l/2+J2^ 



(48) 



which is the proper expression to be compared with Eq. (8). 

As concerns the string equation one comment is in order. In the regime (29) 
the string equation (21) becomes 



(m-n)/2^n 



m/2 



(m - 1/2) 



3 ( ^ 

m-2 —;=i 



u 







(49) 



Recalling that Legendre polynomials satisfy -Pn(l) = 1 we conclude that the solution 
u = y/JI has been consistently chosen. 

2.3 Continuum/discrete dictionary 

Suppose that the m-th critical Hermitian 1-matrix model corresponds to a (p, q)- 
minimal model coupled to gravity. In order to determine p and q we compare scaling 
exponents in the t and regimes. 

In the continuum approach the partition function scales as Z ~ as 
stated in Eq. (10). On the other hand we learned that t should be the coupling 
of the minimal-weight operator, so that in a t-type regime we expect the partition 
function to scale as Z ~ i-Q/Pmin ^ where /3min is the dressing charge of the minimum 
weight operator. These exponents can be written in terms of p and q, 



Q 



P. 



2(p + g) 

{p + q)-l 



a+ p 



(50) 
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Both cases were also studied in the discrete m-th critical model, where the free 
energy scales as JF ~ ^2+1/m jr ^ -j^ ^^le t and regime respectively. 

Identifying the scaling exponents from each approach we find two equations, 

2(p + g) 2 + 1 



{p + q)-l m .g^. 
q 1 

P 2 

which uniquely determine (p, q) = (2, 2m — 1). 

Next we examine the scaling operators. In the (p, q) model there are ^{p — 
l)(g — 1) primary operators, usually labelled Orr', 1 < r < q — 1 and 1 < r' < p—1, 

with conformal weights Arr' = ('"^"'^ '^lpq^'^~^^ ■ With this notation the index v 
that characterizes the wave function {Orr' ) is given by 



. (82) 

For {p, q) = (2, 2m — 1) we have r' = 1 and 1 < r < m — 1 (we can ignore the 
range m < r < 2{m — 1), which only doubles the spectrum of weights). Therefore 
u = m — r — 1/2. On the other hand the discrete model predicts u — n + 1/2. 
We conclude that r = m, — 1 — n, with n = 0, • • • , m — 2, spans the set of primary 
operators and we can finally draw the identification-^ 

?n ^ J d'^e^-'^Om-i-n,! , n = 0, 1, • • • , m - 2 . (53) 

Concerning the remaining scaling operators, we have two special cases: as indicated 
by Eq. (39), am-i corresponds to the boundary operator 

d^ e'*+'^/2 , (54) 

while dm is part of the energy operator, as suggested in Ref. [1] : 

^ (2m -3) ^ 

^rn - J7. r-7T/^o-m-2 ^ energy . (55) 

(zm + 1) 

For the operators an, n > m, we verify that their conformal weights differ from 
the weights of an, n < m, by integers, and thus correspond to secondary operators. 
The 1- and 2-point functions found in the discrete and continuum approaches are in 
perfect agreement. Higher-order correlators involve the issue of fusion rules and lay 
beyond the approximation assumed in Eq. (10). Nevertheless the scaling factors 
agree with the predictions of the continuum. 

Now we are ready to generalize the previous computations for the non-critical 
superstring. 
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3. Supersymmetric non-critical strings 

3.1 Super- Liouville approach to non-critical superstrings 

The results in this section have been obtained in Ref. [2], using the formulation of 
Ref. [11] (see also the review [10]). The total supersymmetric action is given by 
S = SsL + Sm + Sgh, where 

SsL = I E {^-D^^slD'^^sl - ^Y^SL - 4z7Ie"+*--^ (56) 

is the super-Liouville action, representing the supergravity sector of the supcrstring 
in the superconformal gauge; $5^ = + 6'';/' + ^';/' + 99 F is the Liouville superfield; 
/I is a bare cosmological constant; E is the super determinant of the super-zweibein 
and Y the curvature superfield^'^°'^^. Also, 



J E I-D^^mD'^^m - iaoY^Mj (57) 



corresponds to the c < 3/2 supermatter action in the Coulomb gas formulation, 
with matter superfield $m = ^ + + + 99G. The term Sgh stands for the 
ghost action whose explicit form we do not need here. The background charge ao 
is related to the matter central charge: c— f (1 — Sctg)- To compute Q and a+ we 
have to impose respectively that the total central charge ct vanishes and e"+*^^ 
be a (1/2, 1/2) conformal operator. We need the bosonic piece of the super energy 
momentum tensor which is given in components by 

TSL = -1: dct>dct>: -^:,pdiP: +^d^ , Tm ^ -\: dXdX: (dC: -iaod^X 
Tgh = nc + T/3j=:cdb:+:2{dc)b:~:{d^)P:~:^dp , (58) 

where 6, c (/3, 7) are ghost fields. We use free propagators for all fields. The central 
charges for super-Liouville and ghosts are computed in the usual way, that is, csl = 
|(1 + 2Q^) and Cgh = —15. Imposing that csl + c + Cgh = we deduce that 

Q = 2^/^4 + 1. The conformal weight A of an operator e"* is A(e"'^) = -M^iQ) . 

Requiring that A = 1/2 we get a± = —^±^\/Q'^ — 4 = — ^ilctol, with q;_|_Q!_ = 1. 

The particle content consists of a scalar (Neveu-Schwarz, NS) and a spinor (Ra- 
mond, R) vertex, both massless. The Neveu-Schwarz vertex is the supersymmetric 
extension of the tachyon. 



= J d2ze*'=^+^'^ [(z/cC + Pi^) {ikC + 0^) +(3F + ikG] . 



(59) 
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As before, we impose that e*'^*^^^*^^ be a (1/2, 1/2) conformal operator, 

A(e^/c*^+/3*-) = l/c(/c-2ao)-§(/? + Q) = ^ , (60) 

which fixes the NS dressing charge; as in the bosonic case, - ao\- 

The auxihary fields F and G appear in a trivial way in the action. Their 
propagators consist of delta functions, {F{zi)F{zj)) — {G{zi)G{zj)) ~ 5^'^\zi — zj). 

The contractions of other fields typically give (e°"^(^*)e'"^'^^^^) ~ \zi — Zj\ This 
problem will be circumvented by discarding such fields, i.e. we fix F = G = 0. This 
amounts to assuming ab < 0^'^. After such simplification the NS vertex becomes 

ii;^^{k) = ^d^^d^^ e^^^^+^'^«^ = J dz{ikC + pi/j) J dz{ikC + 0i^)e'''''+f^'^ . (61) 

In order to calculate correlation functions we must consider the residual OSP 
(2, 1) symmetry of the superconformal gauge. After choosing zi — 0, Z2 — 1, 
^3 = oo, ^2 = ^3 = and renormalizing the cosmological constant Jc ^ fi, one 
finds^ 

-4N = |^M''+'-'nAQ + l(;3|-fc])) , (62) 

with s = —{Q + '^i/3i)/a+, i.e. the results are completely similar to the bosonic 
case (recall Eq. (8)). 

For the Ramond vertex we have to consider other fermionic contributions. First 
we bosonize the fermions i/j and C into a bosonic massless field h in the usual way, 
that is^ 

ip±iC = V2e^'^ (63) 
with the contraction {h{z)h{w)) = — ln(2; — w) and the superselection rule 

(64) 



(^^e^^'^^"^)^ ^ iff 



Because of supersymmetry the Ramond vertex must also contain a ghost spinor 
field E whose bosonized form T,±i/2 — e^'^^'^ involves a massless bosonic field a with 
propagator {a{z)a{w)) = — ln(2; — w). The two solutions have different conformal 
dimensions, A(e~^/^) = 3/8 and A(e^/^) = -5/8. The requirement 

11 1 

^ Lhh{z)+ikXiz)+m-)^(^z)) = - + -k{k - 2ao) - -PiP + g) + A (S) = 1 (65) 
\ / 8 2 2 

is enough to determine the proper choice. Indeed the Ramond vertex V^ {k, e) should 
represent a massless particle. The on-shell condition {k — cto)^ — {(3 + — 0, 

13 



equivalent to E'^ — = 0, implies that A (S) = 3/8, which selects the solution 
S = S_i/2 = e"'^/^. Therefore our spinor vertex can be written as 

V^{k,e) = V_i{k,e) = yd2ze-^^^^)+^^'^(^)+»'=^(^)+^^(^) . (66) 

From the Dirac equation we find for the dressing the expression^ P{k,e) — — ^ + 
e(fc — tto). There are further versions of these operators obtained by a certain pro- 
cedure defined in Ref. [13], leading to the so-called pictures of vertices. We remark 
that the field a has background charge —2, which implies another superselection 
rule. 



^gg.^(^.)Wo iff = -2 . 

\ i I i 



(67) 

Another useful picture of the NS vertex is the following 

The above vertex is BRST- invariant and can be "picture-changed" into ipNS = 
Qbrstj^'^'n'P ' where ^ is the ghost zero mode (see [13] for details). This is 
lowever not BRST-exact. 

A mixed A'^-point correlator 

/ n N \ 

A^-'^-) = /l[V_,^{k„e,) n i^.sikj)) (69) 

j=n+l I 

is obtained by integrating over matter {Xq) and Liouville ((/)o) bosonic zero modes. 
The total momentum conservation law (7) still holds, supplemented by the rules 
(64) and (67). As in the bosonic case we find a factorizable result^ 

i=l ^ ^ j=n+l ^ ^ 

where the exponent s, defined in (8), now includes spinor vertex momenta. Notice 
that the inclusion of the vertices F_i does not alter the scaling of the partition 

function and Eq. (10) still holds. The insertion of a NS or R vertex modifies the 
scaling behaviour by a factor iJ,~^/°'+ , (3 being the corresponding dressing charge. 

The bosonic case has taught us how important the analysis of 1- and 2-point 
functions are in the comparison between continuum and discrete results. So they 
must be in the supersymmetric theory. As in Eq. (11), momentum conservation 
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and the extra rules (64) and (67) imply that the "area" operator V'jvs(2Q!o) is the 
only vertex with non-vanishing expectation value, scaling as 

(V^,,(2ao))~A^-^/"+-' . (71) 

In the NS sector we have an orthogonality property similar to the one valid before 
for the operators (6), namely 

Concerning the R sector, the superselection rule (64) implies that every 1- 
point function {V_i) must vanish. From Eq. (70) and the rules (64) and (67) 

we also conclude that ^2^''^'' = 0. However, considering the picture-changed NS 
vertex (68) we can build a non-vanishing 2-spinor function, (V_ 1 (fc, e) V_ 1 (2ao — 

/c, — e)';/'^^ ^ (0)) ^ 0. The operator V'^s^(O) = j([^ze~'^'^°^+'^ acts as a screening 
required by the rule (67), its "engineering" dimension is [length]^ and its coupling 
constant must tend to jj? in the //-regime. The screening vertex f^'^ipj^g (0) thus 
defined contributes with a scaling factor /t-"+/"++2 = ^. We therefore have the 
following screened 2-point function 

1 (/c, e) V 1 (2ao - K -e)^i'i;';^' (0)) ^ . (73) 

We could interpret the above result in the following way: the spinor vertex V_i{k, e) 

contributes with the scaling factor while its "conjugate" vertex in this 

representation is the composite operator fA^i/j''^^^ (0) x y_ 1 {2ao — k, — e), which scales 

as This interpretation will be further supported by the discrete model 

results. 

A delicate issue concerning the comparison with the matrix models results is 
the definition of the matter sector of the Ramond vertex, since the gravitational 
and matter fermions merge into e^^, and an identification of matter contribution 
with the Kac table results turns out to be difficult. Recall that in general for a 
Ramond vertex | + | + ^k{k — ao) — \(3{[3 + Q) = 1, where the 3/8 comes from 

the ghost contribution, and 1/8 from e^'*. Therefore, since the latter is neither 
pure matter nor pure gravity, we have to disentangle their contributions and find 
the relation between the pure matter conformal dimension ^Kac {^k) and the 
dressing. Gravitational dressing gives a contribution —^(3{(3 + Q) + jq, as argued 
in Ref. [11]. (The Neveu-Schwarz field is expanded in half- integer components and 
displays no zero mode, while the Ramond field is expanded in integer modes: the 
zero mode contributes 1/16). Put in another way, in order to maintain the form of 
the Virasoro algebra we must have, in the Ramond sector, shift Lq into Lq + j^, 

where 1/16 is the ground-state energy. Therefore we have, Ak — A + j^, where 
A = lk{k - ao). 
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3.2 Discrete super-eigenvalue model 

In Ref. [3] a discrete model was proposed as a supcrsymmetric extension of the 
effective eigenvalue theory, which in its turn comes from the angular integration 
of the Hermitian 1-matrix model. It is defined by the supersymmetric partition 
function 

N N 



n=l i<j 

oo 

VsiK, Oi) = Y,{gkX^, + ^k+i/2diX^i) , 



(74) 



k=0 



where Xi,$i are commuting and anticommuting eigenvalue variables, respectively, 
and gk,^k+i/2 are the commuting and anticommuting coupling constants. 

The planar superloop equations and the double scaling limit of this model 
were studied in Refs. [3,14]. Higher genera were also considered in Ref. [14] but it 
was only after Ref. [15] that a non-perturbative definition was given in terms of a 
supersymmetric KdV hierarchy. 

To compare results with the continuum approach the planar approximation will 
suffice. The m-th critical double scaling limit is similar to the bosonic case. For 
even bosonic potentials the resulting theory is described^ by the supersymmetric 
free energy 

:Fs = -7^dr^[l-dtT+T_dt]u , (75) 

where the functions u and r±, respectively bosonic and fermionic, satisfy the equa- 
tions 

t = U^-J2tnU^ , T±(|.) = ^T±IX" , (76) 

n>0 n>0 

in terms of the renormalized bosonic (fermionic) coupling constants tn (r^)- It is 
also convenient to define the functions Pn[u] = (1 — dfT-^-T-dt)^, generalizing the 
monomials Rn[u] — u^/n of the bosonic theory, whose fiows are given by 

d d d d 

-^Pm[u\ = -^Pn+m[u] , ^-^Pm[w] = Pn+m[u\ ■ (77) 

Otn Ot OTn 0'T± 

Using the string equation (23) we can also write the functions pn as 

p,[«] = (l-D+D_)— , (78a) 

n 

where 

i>0 * 
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In this case we can make explicit the quadratic dependence^'^^ of fermionic couphngs 
exhibited by the free energy: 



rs = {i-D+D_) 



2k 



2^t 



^{1 - D+D_)Tbosonic ■ (79) 



We associate scahng operators to the couphngs, an ^ ^ af~ their 

correlation functions follow easily 



n 



llaa,\ = -^dr^Pa+i[u] , a = ^ai , (80a) 

U=l / i 

^iIl''^) = -^dt-'-^Pk+a+i[u] , (806) 



i=l 



1 d d 



In analogy to the bosonic case and considering the results of Ref. [15] we shall 
call this set of operators/couplings the super- KdV frame. Now we can compute 
correlation functions in the t-regime given by the limits tn,T^ — ^ 0. The free 
energy scales as in the bosonic model, 

jr ^t2+l/m ^ (81) 

while higher-order correlation functions scale as^ 

/ n \ 

n^aO-^'^-^^^-"'^ > (82) 



i=l 

implying the scaling dimensions da„ — ^ and d,,± — ^ Further correlators 

''- lit Tt Tib ^TTl 

vanish. 

We can foresee difficulties with the super-KdV basis similar to the bosonic case: 
indeed we have infinitely many non-vanishing 1-point functions, (ai) ^ 0, while the 
2-point functions {(JiCTj) and {I'fi'J) are not diagonal. The results (71) to (73) show 
that we cannot associate the operators (ui, Vi) to (V'jvs; ) in a simple way. 
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3.3 Superconformal frame of scaling operators 

Let us begin to investigate the dimension of the couphng constants. As in the 
bosonic case the specific heat u has dimension of inverse length, [u] = [length]"^. 
From the string equation and the renormahzation procedure used in Ref. [3] we 
verify that 

[t.] = [length]--"^ , [r±] = [length]--'-^^ . (84) 

Recalling the super-Liouville action (56), wc observe that the "area" operator actu- 
ally has the dimension of a length. In this case the "physical" cosmological constant 
IJ, has dimension of inverse length. We conclude that in general tm-i is the coupling 
with the proper dimension of a cosmological constant rather than t. The correct 
comparison with the continuum results requires that we replace the t-regime by a 
//-like scaling regime. 

As outlined in the bosonic case, we take the string equation (21), consider 
an analytical transformation from the super-KdV basis {tmT^} to a generic one 
{tn, T^}, and replace the t-regime (20) by the following one: 



^m— 1 ^ 

U ^ , 



n ^ m — l,m 



Vn 



(85) 



with solution u = jj, for every m-th critical model. Since we are going to concentrate 
on the properties of 1- and 2-point functions we must make an extra shift of the 
non-vanishing couplings, tm-i — > + tm-i and — > — 1 -|- tm, and then treat 
tn, as perturbative couplings. Imposing also that the dimensions are preserved, 
i.e. \tn] = [tn] and [t^] = [r^], the lowest-order transformations must have the 
form 

oo 

tn = S„//— + , (86) 



with 
and 



Bm-i = -A^^i ; Bn = , n>m (87) 



= E^i"^^V^r±+, . (88) 



s=0 



The critical regime is therefore defined by the limit tn,Tn 0. As we shall see, 
the choice of the same coefficients A^J^^ in both transformations (86) and (88) im- 
plies that (d'iaj) and (pfi^j) are simultaneously diagonalized and also preserves 
supersymmetry. 

From Eqs. (86) and (88) we obtain the new set of scaling operators. 



an = ^A^^^n'an-s 



s=0 



s=0 



(89) 
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In order to find the proper coefficients A^J^^ we sliall follow the orthogonalization 
procedure described in the bosonic model. Prom the transformation (89) and the 
general correlators given in (80) we find 

{^^t^-) = {^i^J) = -^^^'+'+'9ij{A) , (90) 

with the symmetric matrix 

(91) 

Once more we have reduced the issue of diagonalizing 2-point functions to a problem 
of orthogonal polynomials. Notice the differences with respect to the purely bosonic 
case: the integral goes from to 1 now, and the polynomials have no definite parity. 
The solution is given in terms of the Ti-polynomials defined as 



s=0 



(92) 



for which gij = 2^^ij- The corresponding coefficients are given by 

.(m-l) _ Am) .(m) 

^m-s- 2m ' 2m • ^^^^ 

Also the choice (92) guarantees that the only non-vanishing 1-point functions are 

^ , ^^-/-^2m(2m + l) ' ^ ^ 

the former corresponds to the average "area" to be compared with Eq. (71), while 
the latter is expected to represent part of the average energy, as in the bosonic case. 

Therefore we propose the following definition of the superconformal basis of 
scaling operators 

_^ (-l)- {2n-s)l 

h ^' [(n-.)!f^^"^ ' 



T^Q [{n-sy] 



= E 



which is our main result. We emphasize, though, that Eq. (95) corresponds to a 
lowest-order transformation. It can also be inverted, so that the super-KdV frame 
reads 

^ (2n+l-2^) 
= ^ .!(2n+l-.)! ^ ' 

(96) 

r^ = WfY (2n + l-2.) ^ ^ 
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We also notice that, within the approximation (86)- (88), we have 

and the free energy can still be written as = |(1 — ]D^]D-)J^i,osonic- In the new 
regime we verify that the free energy scales as 

while general A'"-point functions scale as 

na,^~//^+^/'+s/^ , (99) 

with the scaling exponents 6^ = n — m and 6 + = n — |, to be compared with 
the continuum values —(3/a-s^. 

The m-th critical string equation becomes 



m—n„,n 

~ 2m 



u \ I u 



= , (100) 



in agreement with the solution u ^ ji due to the general property 7rn(l) = -Pn(l) = 1 
of the TT-polynomials. 

3.4 Super- Liouville/ eigenvalue dictionary 

As in the bosonic case, we suppose that the m-th critical supersymmetric model 
corresponds to some dressed (p, g)-minimal superconformal theory; we determine 
the numbers p and q by analysing the scalings in two different regimes, namely the 
t- and //-averaged regimes. 

In the continuum theory the partition function behaves as Z ~ i-Q/Pmin 
Z ~ . For the supersymmetric (p, q) model, we have ckq = cind 

Q _ 2(P + 9) Q . (101) 



{p + q)-2 a+ p 

In the discrete theory we have found J-g ~ /;2+l/"^ ^^d J-g ~ /z^™"'"^, thus implying 
the following pair of equations: 



2(P + 9) _2+ 1 



(p+q) -2 m 
Q 
P 



1 + - = 2m + 1 



(102) 
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whose solution^ is (p, q) = (2, 4m), m = 1, 2, • • •. The corresponding central charges 
are c 



3 -, _ {2m-lY 
2 ^ m 



0, — • • •. The case m = 1, c = corresponds to pure 



supergravity and is the only unitary model within this series. 

Concerning the spectrum of primary operators, we have two sectors (NS and R): 
the total number of primaries is — l)(g — 1) + 1), which makes up 2m operators 
(m in each sector) in our case. The conformal weight A^^/ of such operators is given 
by the formula 

{rp-r'qf -{q-pf 1 - (-1)^-^' 
^^'" = 8^^^ + 32 ' 

with 1 < r' < p — 1 (i.e. r' = 1 in our case) and l<r<g — 1 = 4m — 1. When the 
difference r — r' is even, the operator belongs to the NS sector. 



NS i{i — 2rn + 1) 

4m 



^2i+i,i = ^ — ' i = 0,l,---,2m-l , (104) 



and for r — r' odd one has the R sector 

(i + 1 -m)2 - (m- 1/2)2 ^ 



^^+2,l= ' ^ ""'4^'"" +16 ' * = 0,l,---,2m-2 . (105) 



In fact we can restrict the label i to the range i = 0, • • • , m — 1, since the remaining 
values only duplicate the spectrum. Prom the weights A^r' it is straightforward 
to calculate the corresponding momenta krr> and dressing charges (3rr'- Then, by 
comparing the values of —Prr' / Oi+ with the exponents in Eq. (99) we can associate 
the bosonic scaling operators with the NS vertices, 

^ ^Ns{K(m-n)-i,i) ' n^O,---,m-l , (106a) 
and the fermionic operators with the R vertices in the following subtle manner 

^ V^-i(W.).i) ' (1066) 
^ /.2vL?(0)Vli(2ao-A;,(_,,J , n = 0,---,m-l . (106c) 

The special operator dm-i corresponds to the "area" or cosmological term. We can 
also define an energy operator, 

(2m -1) ^ 

Sm = (rm- (2to+1) ^'^'"~^ ' ^ ^ 

which points out the role of the operator a^. The remaining operators correspond 
to secondary fields. Notice the absence of a boundary operator, as opposed to the 
bosonic case. 
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3.5 Wave functions and minisuperspace approximation 

We shall determine the supersymmetric version of Eq. (15) obeyed by the wave 
functions of scaling operators in the superconformal frame. In this way we expect to 
circumvent some ambiguities, which usually arise in attempts^^ to derive the super 
WdW equation in minisuperspace approximations from the continuum theory. 

We start from the macroscopic superloop W±{1, 6±) defined in Ref. [3], whose 
expectation value reads 

{w±{i,e±)) = ^{-di' + T+T_dt±e±d^\^dt)e-'^ . (108) 



Using the string equation and the definition (786), this equation can be rewritten 



1 

{W±{l,e±)) = [l-D+D_^e±D^]^= dye-^-(^) , (109) 

AvVTTi J to 

to be compared with Eq. (23). In fact we have the same integral of the bosonic 
model, implying an expansion in terms of the functions i(^n defined in Eq. (25). 
Using the following identity between Bessel functions and 7r-polynomials, 

e-/^^„«4i)-y!fd-„©e-' , (110) 

we can calculate the integral in Eq. (109) in the superconformal frame, finding as 
a result the following expansion 



kx/ttI Jt( 



oo 1 m+1/2 



+ ^f;f„«-+V2e-W2^^^^/JM . (ni) 
We therefore obtain, in the regime tn,T^ 0, the following wave functions 

(112) 



KTT \ ^ J 

Both satisfy the wave equation 



4 + aV^i!i! + ('i + iVi*. = o, (113) 



dl 2 ; 4 V 2 
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which we expect to be the supcrsymmetric WdW equation in the minisuperspace 
approximation. Indeed, after dropping the space dependence of the fields in the 
continuum theory, we have the super-Liouville Hamiltonian 



suggesting^^ the associations 

p2 + ^^^e^-^<^ ^ _ (^/| + , (115a) 



A'' „a+<b , ..2;2 



2a% 



^ . (1156) 



Notice that the fermionic contribution is summarized by an /-dependent shift of the 
momentum p, which amounts to the extra exponential factor e~^'^/^ in comparison 
to the Bessel wave functions of the bosonic model. The understanding of this effect 
in terms of the continuum formulation is at present under investigation. 

4. Conclusion 

Even though a matrix formulation is still missing, we have found the superconfor- 
mal background of couplings and operators of the discrete super-eigenvalue model, 
which can be completely translated into the super-Liouville language, and a pre- 
cise correspondence to the dressed (2,4m) minimal superconformal theory can be 
established. As a by-product we have derived the wave equation that should be 
equivalent to the proper minisuperspace approximation of the WdW equation that 
supcrsymmetric loop amplitudes are expected to obey. 
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